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ABSTRACT 

We present two new types of measurements that can be made with diffusing-wave spectroscopy, a form of 
dynamic light scattering that applies in the limit of strong multiple scattering. The first application is to measure 
the frequency-dependent linear viscoelastic moduli of complex fluids using light scattering. This is acc omplishe d by 
measuring the mean square displacement of probe particles using DWS. Their response to thermal fluctuations is 
determined by the fluctuation-dissipation relation, and is controlled by the response of the surrounding complex fluid. 
This response can be described in terms of a memory function, which is directly related to the complex elastic 
modulus of the system. Thus by measuring the mean square displacement, we are able to determine the frequency 
dependent modulus. The second application is the measurement of shape fluctuations of scattering particles. This is 
achieved by generalizing the theory for DWS to incorporate the effects of amplitude fluctuations in the scattering 
intensity of the particles. We apply this new method to study the thermally induced fluctuations in the shape of 
spherical emulsion droplets whose geometry is controlled by surface tension. 


A microgravity environment has significant potential benefits for a wide variety of scientific problems 
involving colloidal particles or other complex fluids, where sedimentation can play an important role. Several 
experiments are currently planned to exploit this environment Laser light scattering is an important experimental 
technique for studying the physics of these complex fluids; as a result a laser light scattering apparatus is being 
constructed by NASA to allow measurements of different types to be made in space. This apparatus will allow the 
traditional light scattering techniques to be employed, including both static and dynamic light scattering, hi addition, 
the apparatus will allow several ancillary measurements to be performed. Among the most important is the 
measurement of the rheological properties of the complex fluids being studied. Currently this measurements will be 
done by mechanically agitating die sample cell and measuring the response of the complex fluid. This will provide a 
measure of the complex, frequency-dependent linear viscoelastic modulus of the suspension, G*(oa) - 
G'(o))+iG"((o). Here, the real part, G'(ffl), reflects the energy storage or the elastic component of the modulus, 
while the imaginary part, G"(ta), reflects the viscous dissipation of the modulus. While these mechanical 
measurements of the modulus provide a direct measure of the viscoelastic response of the fluid, they are limited in 
the frequency range that is accessible. Moreover, since the primary function of the NASA apparatus is light 
scattering, a preferable method for measuring the visco-elastic response would be to use light scattering techniques. 

In this paper, we discuss a new method for measuring the viscoelastic response of a complex fluid using 
solely light scattering methods. We discuss the use of diffusing-wave spectroscopy (DWS), which is form of 
dynamic light scattering applicable in the strongly multiple scattering limit[l,2]. Two new applications of DWS are 
presented. We first discuss the use of DWS, or more traditional dynamic light scattering techniques, to measure the 
frequency-dependent linear viscoelastic modulus of complex fluids. We also discuss the extension of DWS to include 
amplitude fluctuation due to the fluctuating shape of the scattering particles, and use this to measure the shape 
fluctuations of emulsion droplets, where surface tension controls the spherical shape of the particles. 

We begin by presenting a novel method for measuring the linear viscoelastic properties of a complex fluid 
over an extended range of frequencies. We show that the response of the fluid to thermal fluctuations, as probed by 
the average motion of small particles dispersed within the fluid, provides a very close representation of the response 
of the bulk fluid to an imposed shear strain. The essential physics of this approach is that the bulk mechanical 
susceptibility of the fluid appropriately describes the microscopic viscoelastic forces acting on a small particle 
excited by the thermal stochastic forces which lead to Brownian motion. By describing the stochastic motion of the 
small probe particle with a Langevin equation generalized to include viscoelasticity and imposing energy 
equipartition using the fluctuation-dissipation theorem, the average particle motion, given by the mean square 

displacement, (Ar 2 (t^j, determines the time-dependent memory function, £(0, of both the surrounding and bulk 

complex fluid. This time-dependent memory function is proportional to the more commonly encountered stress 
relaxation modulus and contains the same information as G*(cu). We use dynamic light scattering to measure 

(Ar 2 (tfj , and compare the G*(co) obtained with that measured directly by mechanical means using a rheometer. 

We describe the motion of a small, neutrally buoyant particle dispersed in a complex fluid by means of a 
generalized Langevin equation, 
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mv(t) = fait) - - T)v(r)dT 


( 1 ) 


where m is the particle mass and v(r) is the particle velocity. This equation expresses the balance of forces on the 
particle, relating the stochastic Brownian forces, fait), to the response of the surrounding fluid, which is described by 
a generalized time-dependent memory function, £(f), and reflects the dependence of the properties of the material at 
time t on the behavior at previous times. A simplified form of this equation is commonly used to describe the 
motion of a particle in a purely dissipative viscous fluid; here we generalize it by using the memory function to 
describe the average, effective-medium response of a complex fluid that can also store energy. Allowing the medium 
to store energy profoundly changes the temporal correlation of stochastic forces acting on the particle at thermal 
equilibrium. Because energy from thermal fluctuations can be stored in the medium over time scales determined by 
the memory function, the dissipation required to maintain thermal equilibrium also occurs over these time scales. 
As a result, the stochastic forces acting on the particle are correlated over finite rather than infinitesimal time scales 
to prevent an unphysical energy buildup. Thus, the fluctuation-dissipation theorem differs from the commonly 
encountered delta-function correlation of a purely viscous fluid, and here takes the form[3]: 

(/b(0)/b(0) = WW (2) 

where the coefficient ensures thermal equilibrium, with kg Boltzmann's constant and T the temperature. 

By taking the Laplace transform of Eq. (1), and using Eq. (2), the viscoelastic memory function can be 
related to the velocity autocorrelation function, and, hence, the mean square displacement of the particle. We assume 
that the Laplace transform of the microscopic memory function is proportional to the bulk frequency-dependent 

viscosity, ij(s)= £(s)/6jca, where s represents the frequency in the Laplace domain. We choose the coefficient to 


ensure the correct Stokes drag on a particle of radius, a, in a purely viscous fluid, recognizing that this coefficient 
may differ somewhat in describing the behavior of a particle in a viscoelastic fluid. The frequency-dependent shear 

modulus is directly related to the viscosity by G(s) = srj(s ) , and hence to the mean square displacement, r 2 (r)^: 






-ms 


(3) 


The terms within the brackets come from the solution of the Langevin equation for the memory function, £($); the 
first captures the thermal fluctuation-dissipation in the medium surrounding the particle, while the second is due to 
the particle's inertia. In most complex fluids, the inertial term is negligible compared to the fluctuation-dissipation 

term, except at very high frequencies. Fitting the optically measured G(s) with an appropriate functional form, die 
complex, frequency-dependent shear modulus can be determined using analytic continuation by setting s = ico, and 
identifying G'(<a) and G"(«o) as the real and imaginary parts. Provided the functional form describes G(s) ova* all 
measured frequencies, the Kramers-Kronig relationship will be satisfied over the measured range of frequencies. The 
complete procedure establishes a very general relationship between the mean square displacement of suspended 
particles and the bulk rheological properties of the complex fluid surrounding them. 

To test the applicability of this scheme, we use diffusing-wave spectroscopy to determine the mean square 
displacement of the particles in a concentrated suspension of silica particles of relatively uniform radius, a = 0.21 
/rm, suspended in ethylene glycol. These particles behave essentially as hard spheres. Their volume fraction was 
increased by centrifugation to <p~ 0.56, quenching in disorder, to form an isotropic hard sphere glass. Diffusing- 
wave spectroscopy was used in the transmission geometry to measure the correlation function for a 4-mm thick 
sample, shown in Fig. 1; characteristic of a colloidal hard sphere glass, it exhibits an initial, rapid decay to a plateau 
value, followed by a final decay at longer times. The particles are relatively small, so that DWS is slightly sensitive 
to collective motion; nevertheless to a good approximation the correlation function can be inverted to obtain the 
mean square displacement of the particles. This requires knowledge of the transport mean free path of the light, /*, 
which was obtained from static transmission measurements. The mean square displacement obtained from the 
correlation function is shown in Ftg. 2. The linear increase at the shortest times reflects the initial diffusive motion 
of the particles; at longer times their motion is constrained by the local cage structure of the colloidal glass. The 
decay in the correlation function at the longest times can not be simply interpreted as a mean square displacement of 
individual particles setting the upper bound on the data accessible. 

We numerically calculate the Laplace transform of the mean square displacement, and use Eq. (3) to 
determine G(s), which is shown by the open points in Fig. 3. The accuracy of the Laplace transform is limited at 
the shortest and the longest time scales because of the limits of the available data. Nevertheless, the data exhibit 
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several very general features that characterize G(s). The region having relatively low slope at lower frequencies 
indicates energy storage and hence a range where the real component of the complex modulus dominates. The region 
having nearly unity slope at higher frequencies indicates viscous dissipation, and hence a range where the imaginary 

component of the complex modulus dominates. This is very generally true of the data for G(s); regions of small 
slope indicate that the elastic component of the modulus dominates, while regions of near-unity slope indicate that 
the viscous component dominates; the deviations from these slopes reflect contributions from both components. 

To obtain the real and imaginary parts of the complex modulus, the frequency-dependent magnitude given 

by G(ico)G(-ico) = |g*(£ 0)| can be fit using a superposition of complex moduli which satisfy the Kramers-Kronig 
relationship and have a single Maxwellian relaxation time. Rather than follow this more general procedure, we 
instead use our physical intuition about the behavior expected for a colloidal hard sphere glass to determine G(s). 
We include a constant to account for the elasticity at low frequencies; a term proportional to 5®-^ to account for the 
behavior of the plateau, as predicted by mode coupling theories[4]; a term proportional to to account for the 
predicted asymptotic high frequency elastic modulus of hard spheres [5]; and a term proportional to s to account for 
the high frequency viscous component. A fit to a sum of these terms yields the solid line in Fig. 3. The real and 
imaginary components of G*(a>) can then be calculated directly. Due to the form of G(s), these do not satisfy 
Kramers-Kronig at low frequencies, but they do properly predict the moduli within the limits of the measured data. 

To compare the elastic moduli calculated from the light scattering data with those determined directly by an 
oscillatory mechanical measurement, we use a controlled strain rheometer with a double-wall Couette sample cell 
geometry. A small, sinusoidal strain was applied and the in and out of phase components of the resultant stress were 
determined as the frequency was varied. The amplitude of the strain was maintained well below 1% to ensure that the 
data were obtained in the linear regime; this was verified by varying the strain amplitude. The upper frequency was 
limited by the motor's ability to produce a reliable sinusoidal strain, while the lower frequency was limited by the 
sensitivity of the torque transducer. The resultant real and imaginary components of the complex modulus are shown 
in Fig. 4. The dashed lines show the values of G'(co) and G"(co) obtained from the light scattering data. 
Remarkably good agreement is found, except for the loss modulus at the lowest frequencies. The light scattering 
data represent the overall trend of the rheological data quite well, particularly for the larger storage modulus. 

These results clearly illustrate the possibility of using light scattering methods to determine the visco- 
elastic response of a suspension of hard spheres. We have also successfully tested the technique on very different 
forms of complex fluids, including polymer solutions and concentrated emulsions; in all cases tried the results are in 
good agreement with the mechanical measurements. The essential physics of this technique relies on the equivalence 
of the frequency-dependence of the microscopic memory function which describes the response of a probe particle in 
the complex fluid and the macroscopic complex viscoelasticity of the complex fluid. Provided particle inertia can be 
ignored and the viscoelastic coupling between neighboring particles varies as (p, the complete equations of motion 
for a lattice of interacting particles can be reduced on average to a single particle equation of motion like the 
generalized Langevin equation having a memory function which describes the average bulk dynamics inherent in the 
microscopic coupling. This equivalence is perhaps somewhat surprising; the motion of a probe particle, as probed 
by light scattering normally reflects relatively large wavevector, q, behavior. Moreover, light scattering typically 
probes the longitudinal response of a system. By contrast, the elastic modulus of a material reflects the transverse 
response of the material in the limit of q -» 0. However, the mean square displacement does couple to shear modes; 
this is certainly true in the limit of a simple viscous fluid, where the shear viscosity determines the particle's 
motion. Thus, by analogy, we expect the shear elasticity of a complex fluid to determine the mean square 
displacement of a probe particle. However, the elasticity is typically measured by applying a uniform strain across 
the material; this is not the case for a small probe particle. As a result, the flow pattern of the viscoelastic fluid 
around the particle may differ significantly from the case of a simple viscous fluid; this will modify the coefficient 
relating the viscosity to the memory function, and may even make it frequency dependent Thus, this method may 
not provide quantitatively exact measures of the elastic moduli; nevertheless, the overall trends are captured, and, as 
shown by our results, the agreement is surprisingly good. 

Perhaps most importantly, this method provides a measure of the elastic moduli over a much broader range 
of frequencies than is obtained with mechanical measurements. Moreover, it provides a convenient method to relate 
the energy storage and loss due to strain to the microscopic mean square displacement of the particles in the 
suspension. These two advantages can provide considerable new insight into the underlying physics of the complex 

fluid. For example, for the hard sphere suspension, the plateau in (Ar 2 (tfj is known to reflect the quenched disorder 

of a colloidal glass; the particle can diffuse freely at very short length scales until it reaches the cage formed by the 
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neighboring particles. Then the particle motion is restricted until the cage decays at much longer time scales. This 
behavior is accounted for theoretically by mode coupling theory (MCI) [4,6]. Our results indicate that the cage 
effects also play an important role in the rheological behavior of a hard sphere suspension; the existence of the cage 
leads to an additional mechanism for energy storage; presumably the strain distorts the shape of the cage, allowing 
energy to be stored, and causing the real part of the complex modulus to dominate. Mode coupling theory accounts 
for this behavior at high frequencies and predicts a /3-relaxation contribution proportional to j®** in addition to a 
constant nonergodicity parameter which characterizes the low frequency plateau elasticity. This component 
improves our fit to the transformed light scattering data, and directly contributes to the frequency dependence of the 
storage modulus. In fact, the predictions of MCT can be used to better describe the behavior of the loss modulus at 
lower frequencies. The rise in G”(co) measured mechanically suggests an additional decay process at lower 
frequencies; MCT predicts that this a-relaxation process, or von Schweidler decay, results from the ultimate cage 
breakup at long times, and provides a functional form to describe it, suggesting an additional term proportional to 
j-0.55 }„ G(s). Subtracting this a-decay term leads to much better agreement at low frequencies between G"(a>) 
obtained from light scattering and that measured theologically, as shown by the solid line in Fig. 4. 

The second new application of DWS discussed here is the generalization of the theory for the technique to 
include amplitude fluctuations of the scattering particles. We use this scattering method to study the thermally 
driven shape fluctuations of monodisperse emulsion droplets, and to determine the effects of increasing volume 
fraction on the fluctuations. The thermally induced increase in the surface area of a fluid droplet is A « IcbT/T, 


where r is the surface tension. For a typical surface tension of 10 dynes/cm, this corresponds to an increased area of 
only about 30 A^; for a 1 pm diameter droplet, this represents a change in the radius of less than 0.1%. The 
relaxation time, r, of these fluctuation is also very rapid; it is determined by the viscosity of the two fluids. If the 
viscosity of the droplet is much greater than that of the continuous phase, the relaxation time is for 

ij = 0.1 P and R » 1 pm, r « 10 ' 6 sec. The combination of the very short time scale and the very small 
amplitude of the thermal fluctuations of emulsion droplets have, to date, precluded their observation. 

Our samples are emulsions of silicone oil in water, stabilized with sodium dodecylsulphate, and purified 
using the technique of fractionated crystallization to yield highly monodisperse droplets with a radius of R » 1.4 
/an[7]. The surface tension of the interfaces is measured to be 9.8 dynes/cm. We perform our DWS measurements 
in the transmission geometry using an expanded beam from an Ar + laser, which gives &o = 16.3 pm"* for the 
incident wavevector in water. The detected light is collected from a point on the exit side of the 5-mm-thick sample 
cell. The transport mean free path, /*, is determined independently by a static transmission measurement [2], 
allowing the correlation function to be inverted to determine the dynamics of the individual emulsion droplets [ 8 ]. 

We might expect the dynamics of emulsion droplets to be the same as those of hard spheres. This is indeed 
the case for an emulsion made of a high viscosity oil, i\ = 1000 cP. This is illustrated in Fig. 5, where the open 

circles represent (a r 2 (r)\, the mean square displacement of an individual droplet The data follow the shape 


predicted for hard spheres, shown by the solid line. By contrast an emulsion made from a much lower viscosity oil, 
7 j = 12 cP, exhibits distinct additional dynamics, as shown for </> = 0.35 by the solid circles in Fig. 5. At early 
times, the data are significantly higher than those expected for solid spheres; at later times the data merge. Similar 
behavior is observed for all <j>. These additional dynamics result from the shape fluctuations; for the higher tj, the 
relaxation time is so long that any shape fluctuations are masked by the large displacements that occur during r. 

To quantitatively describe these new dynamics, we must generalize the formalism for DWS to include 
fluctuations in the scattering amplitude, in addition to the phase fluctuations which result from translational motion. 
The analysis of the correlation function measured by DWS entails the calculation of the contribution of diffusive 
light paths comprised of a large number of scattering events; these are all represented by an angle-averaged scattering 
event [2]. Thus, to include the amplitude fluctuations, we write the correlation function for this ^-averaged, 
scattering event as 




(b*(q,t)b{q,Q)y 

(b*(q,0)b(q,0)) 


U ex p{-iqAr(t)}) . 


(4) 


Here, the time dependent scattering amplitude is b(q,t), and the brackets with the subscript q denote an ensemble and 
g-average. We assume that the amplitude fluctuations are independent of the translational motion of the droplets, and 
we neglect the effects of correlations between the particles [9]. The amplitude fluctuations contribute an additional 
correlation function which is assumed to consist of a sum of a constant portion and a much smaller fluctuating 
portion. Performing the 4 -average, 
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( 5 ) 




where 00 is the total cross section and A cKf) is the correlation function of the fluctuating portion of b(q,t), ensemble 
and ^-averaged. We have also used the relationship, =k$ll 21*, where l is the scattering mean free path [10]. 


To calculate the correlation function for a diffusive path of length s, consisting of n - sll scattering events, 
we take the product of Eq. (S) n times [2]. Since Ao (0« o 0, the first term can be approximated as an exponential, 
exp { (Ao(O-Acj(O))/ 0O } ; this ensures that the contribution of a single diffusive path is still linear in s, and allows the 
standard DWS analysis to be retained. Thus, the total DWS correlation function is a sum over the contributions of 
all paths, weighted by the probability that a photon follows the path, P(s): 


g 1 (0 = J/ , (5)exp| 


ski 

3 1* 


(A r 2 (0) 


3 / Acr(O) - Aa(t) 

Ik l <*o 




( 6 ) 


Since P(s) is known [2], this equation can be solved in exactly the same manner as done for DWS from solid 
particles. The effects of the amplitude fluctuations are contained in the second term in the square brackets. It is clear 


from Eq. (6) that DWS probes the relative fluctuations of the cross section; like 



their contribution increases 


from zero at (=0 but saturates at long times. It is also clear why very small fluctuations can be detected; the signal 
arises from the sum of n independent amplitude fluctuations. These minute fluctuations would not be detectable 
without the advantage of the multiple scattering. 

We note that this derivation is independent of the nature of the amplitude fluctuations; not only shape 
fluctuations, but other phenomena, such as rotational motion of aspherical particles can also result in a similar 
contribution. Moreover, the DWS data can be inverted, and the translational and amplitude contributions can be 
separated. We do this here by subtracting the translational motion; the full effects of the hydrodynamic interactions 
are included by using the scaling form which describes the data for solid spheres [8]. We plot B?Acr(t)la 0 for <j> - 
0.35 in Fig. 6; we multiply the normalized fluctuating cross section by to give the correlation function of the 
fluctuations of the droplet radius. We note that the data exhibit a distinctly non-exponential decay. 

Similar behavior is found for all other droplet volume fractions. Surprisingly, the shape of the correlation 
function is independent of <fr, all the data can be scaled onto a single master curve. This enables us to determine the 
^-dependence of the shape fluctuations. We find that the characteristic frequency of the decay decreases approximately 
linearly with increasing volume fraction, as shown by the open triangles in Fig. 6. By contrast, we find that the 
amplitude of the shape fluctuations increase dramatically with <jr, this is shown by the solid circles in Fig. 6. 

To quantitatively describe the correlation function of die fluctuating shape, we expand the instantaneous 
radius of the drop in spherical harmonics [1 1]. Each independent deformation mode creates an excess area given by 
ksTHT. Together with the conservation of volume, this condition determines the amplitude of the expansion 
coefficients. Each mode relaxes independendy with an exponential decay rate, ft)/, that must be calculated in the 
overdamped, or low Reynolds number, limit [12]. We obtain 


Act (t) = 


k B r V 
4n TR 2 jr 2 


21 + 1 
/(/ + l)-2 


■%i exp{-cu/t), 


(7) 


with 


1(1+2X21+1) r 
1 2(2?+ 41+ 3) nR' 


( 8 ) 


In general, «a/ depends on the viscosities of both fluids [13]; here we have assumed that the viscosity of the oil is 
much greater than that of the water. Note that asymptotically mj-l, reflecting the fact that the dynamics are driven 
by surface tension; by contrast, mi~fi when the dynamics are driven by the interfacial rigidity [11]. The cross 
section for each mode is calculated within the Rayleigh-Gans (RG) approximation [14], 

81 = 7rjfin0(l + cos 2 #)[3 ji(x)f d9 , (9) 


where x - 2&of?sin(0/2), and where we have omitted an optical constant which is canceled in the normalization by 
ctq. The coefficients, g/, depend only on koR, and become negligible when hkgR, since the length scale of the 
features described by such modes is much smaller than the wavelength of the radiation. This sets the upper bound to 
the expansion; we find that the series converges for / » 20. Because of their large mismatch in index of refraction, 
Mie theory, suitably generalized to treat non-spherical scatterers, should be used. However, the average cross section 
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calculated within the RG approximation differs by less than 0.5% from that calculated using Mie theory. Thus, we 
conclude that the RG approximation is sufficiently accurate. 

These expressions apply only to an isolated droplet However, since the functional form of our data is 
independent of </>, we can use them to describe all the data by allowing both the characteristic frequency and the 
amplitude to vary. The solid line through the data in Fig. 6 illustrates the fit for <p =0.35. Excellent agreement is 
obtained. The non-exponential decay clearly reflects the contribution of the different modes. Quantitative agreement 
between the theory and the experiment should be obtained when the ^-dependent characteristic frequency, <a($), and 
amplitude, Ao(0), are extrapolated to zero volume fraction. Our extrapolated value of <n(0) is identical to the 
predicted value of T/rjR = 1.2 MHz. Similarly, the extrapolated value of the amplitude ratio R^Aa(Q)/a 0 is 65 A 2 , 
in good agreement with the predicted value of (*fir/4M'Oo)2(2/+l)g//(/(/+l>2) = 52 A 2 . The remaining discrepancy 
may reflect the error introduced by the RG approximation. 

The ^-dependence of both the frequency and the amplitude must reflect the consequences of the interactions 
between the droplets, and these are not, as yet, well understood theoretically. The normalized characteristic frequency 
is well represented by a linear behavior, ra(0)/ra(O)=l-O.780. The frequency depends only on the surface tension, R, 
and the viscosities of the two liquids; it is unlikely that T changes with <j>, so that the most likely origin of this 
behavior is the ^-dependence of the effective viscosity of the emulsion. A similar linear form was recently predicted 
for the 1=2 modes, but the predicted coefficient is 1.4 [15]. However, this theory was restricted to the case of equal 
viscosities of the oil and water, which may account for the discrepancy. 

The pronounced ^-dependence of the amplitude of the shape fluctuations is more surprising. This may 
reflect the effects of collisions between the droplets, which become increasingly likely as <f> increases; because the 
droplets are flexible, these collisions may result in additional deformations. We model the effects of collisions by 
introducing a contribution to the amplitude proportional to the fraction of the droplets that, at any instant, are 
colliding with a neighbor. This fraction we estimate as 0g(2,0), where g(2,0) is the pair correlation function at 
contact. The dashed line through the data in Fig. 6 is a fit of Ao(0)/Ao(O) to die functional 1 form l+C</>g(2,<j>) with 
g(2,0)=(l-<p/2)/(l-0)3 as holds for hard spheres [16]. Good agreement is found with C=1.7, suggesting that colliding 
droplets deform an extra amount that is about 1.7 times the deformation exhibited by an isolated droplet More 
generally, it should be possible to regard these increased shape fluctuations as the consequence of the osmotic 
pressure of flexible spheres. The expression 1+C0g(2,0), with C=4 and g( 2,0) as given above, is in fact the ratio of 
the full osmotic pressure of a hard-sphere suspension to the kinetic part nky$T, where n is the number density [16]. 
This observation suggests that Ao(0)/Ao(O) may be related to the reduced osmotic pressure of the flexible droplets. 

The results presented here are two extensions of DWS. Amplitude fluctuations can arise from many 
phenomena in addition to the shape fluctuations presented here; for example, rotational diffusion of aspherical 
particles will also lead to similar effects, and the theory presented here should be applicable. The ability to measure 
the linear viscoelastic modulus optically enables its determination over a much broader frequency range, without 
mechanical motion. This could be particularly beneficial for experiments studying complex fluids in microgravity. 

We gratefully acknowledge very useful discussions with T. Lubensky, R. Klein, W. Russel, W. Graessley, 
J. Bibette, S. Milner and W. Cai. 
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FIGURES 




Fig. 1 (LEFT) DWS autocorrelation function for a suspension of hard sphere colloids at 0 •« 0.56. 

Fig. 2 (RIGHT) Mean square displacement of hard spheres determined from the data in Fig. 1. At early times 
die particles diffuse; at later times the particles are trapped in cages formed by the neighboring particles. 




Fig. 3 (LEFT) Laplace transform of the complex modulus, G(s), determined from the mean square 
displacement, shown in Fig. 2, for hard sphere colloids. 

Fig. 4 (RIGHT) Comparison of the complex moduli for hard sphere colloids determined with light scattering 
(solid lines) and with a rheometer (points). The dashed line does not include a decay at longer time; the 

solid line does. The light scattering data are determined from the G(s) shown in Fig. 3. 




*(sec) tQise c) 

Fig. 6. (LEFT) Inverted correlation functions of monodisperse emulsions. The data are the sum of the mean 
square displacement and a term related to the correlation function of the fluctuating scattering amplitude. 
The solid points are for an emulsion with an oil viscosity of 12 cP, while the open points are for an 
emulsion with an oil viscosity of 1000 cP. The solid line is the theoretical prediction for rigid spheres. 
Fig. 6. Correlation function of the radius fluctuations of an emulsion with <p =0.35, compared with theory. 


359 





Fig. 7. 



0 


The ^-dependence of the characteristic frequency (open triangles) and the fluctuating amplitu d e (solid 
circles). Both quantities have been reduced by their extrapolations to zero volume fraction. The solid line 
is a fit to a linear decay, with a coefficient of 0.78; the dashed line is a fit to a simple model which 
attempts to include the effects of interactions between the droplets. 
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